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0> . Abstract 

We review recent numerical progress in the study of finite dimensional strongly 
'^ ■ disordered magnetic systems like spin glasses and random field systems. In particular 

Q I we report in some details results for the critical properties and the non-equilibrium 



dynamics of Ising spin glasses. Furthermore we present an overview over recent 
investigations on the random field Ising model and finally of quantum spin glasses. 
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1 Introduction 

Spin glasses and random field systems are magnetic materials in which a structural 
disorder occurs as a consequence of a special preparation process. The latter either 
changes the chemical composition of compounds and alloys (via dilution or mixing) 
or it decrystallizes the pure material (via sputtering or milling). The effect is in 
any case a randomness in the position of the spins and/or the sign and strength 
of the interactions among them. If this disorder is strong enough it can drastically 
change the magnetic properties of these material and can give rise to a new sort of 
phase transitions. In many cases even a completely new kind of low-temperature 
phase, the spin glass phase, might occur, which is in some sense ordered but does 
not possess a translationally invariant magnetic pattern. Theoretical models for 
these phenomena exist for over 20 years, however, it took about half of this time 
to establish the most fundamental facts for these models like the mere existence of 
a phase transition. The characterization of the latter transitions via the determi- 
nation of the critical exponents is still far from settled. The investigation of the 
low temperature phase and its static and dynamic properties is a matter of ongoing 
experimental and theoretical research. A major role in these activities is taken by 
Monte Carlo simulations, which have found wide-spread applications in condensed 
matter physics |^ . 



Due to the enormous research activity on spin glasses many excellent reviews on 
this topic have already appeared during the last decade. To begin with. Binder and 
Young 1^ and Chowdhury |^ give a complete survey on the spin glass literature 
until the year 1985. Also the proceedings of the two Heidelberg Colloquia on Spin 
Glasses |p.l8| , |119|| contain useful references. Mezard et al. ||181|| focus on the mean- 



field theory of spin glasses and the book by Fischer and Hertz |^ is a rather 



complete introduction into this complex field. Mydosh's book ||190|| presents the most 
recent overview over experiments on spin glasses and a very stimulating collection of 
articles on random magnets ||243|| gives insight into ongoing experimental activities 



on spin glasses and random field systems. Excellent reviews of the latter have been 
given by Nattermann and Villain ||192|| focusing on theoretical aspects, by Belanger 



and Young ||T^ comparing theoretical predictions for the critical properties with the 
experimental situation and by Kleemann ||152|| about experimental evidence of the 
so called domain states, which is a characteristic feature of these systems. 

Hence we do not intend to repeat the effort of the above reviews — we shall 
focus on frustrated Ising systems and its investigation via Monte Carlo simulations 
during recent years. In section 2 we review the status quo concerning the critical 
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behavior of finite dimensional Ising spin glasses. Section 3 is devoted to the at 
present rapidly evolving field of the non-equilibrium, low-temperature dynamics of 
Ising spin glasses and other models exhibiting glassy features. Section 4 surveys 
our present knowledge of the critical properties and the off-equilibrium scenario 
of the random-field model and the diluted antiferromagnet in a uniform field. In 
section 5 we review the numerical methods that are useful in Monte Carlo studies 
of slowly relaxing systems like those considered here. Section 6 discusses quantum 
effects in spin glasses and especially the zero temperature quantum phase transition 
occuring in transverse-field Ising spin glasses. Finally section 7 lists other models 
like Heisenberg, Potts, orientational, vortex and Bose glasses. 

2 Ising spin glasses: Critical properties 

The simplest model that incorporates the necessary ingrediences for a spin glass 



as quenched disorder and frustration is the Edwards- Anderson model pS], which is 
defined by the Hamiltonian (or energy function) 

H=-J:J^AS,-hJ:S^. (1) 

The spin variables are of Ising type, i.e. Si = ±1, which are placed on a (i-dimensional 
hyper-cubic lattice (simple cubic in three dimensions). The interactions among them 
are short ranged, which is indicated by the sum over only nearest neighbor pairs 
{i,j)- The interaction strengths Jij are quenched random variables that obey a given 
probability distribution, most commonly a Gaussian with mean zero and variance 
one or a binary distribution Jij = ±1 with probability one half. The last sum in (|l]) 
describes the effect of an external magnetic field ii h y^ 0. A typical experimental 
realization of this model is the short range Ising spin glass Feo.sMuo.sTiOa ||136| , |114| 



mixing ferromagnetically and antiferromagnetically interacting components. 

A complete survey on the results for the critical properties of the two-, three- 



and four-dimensional Ising spin glass known until 1986 can be found in ||2^. In 
essence the spin glass transition (at temperature Tc) in the infinite system is charac- 
terized by a divergence of the non-linear (and not the linear, as in e.g. ferromagnets) 
susceptibility 
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which is caused by a divergence of the correlation length ,^ ~ (T — Tc)~^ that sets 



{T-T,)-\ (2) 



the characteristic length scale for the decay of spin correlations defined via 

G(r) = [{SA+r)\. - r-^"-'^^^ ~9ir/0 ■ (3) 

Here (■ ■ ■) means a thermal average and [■ ■ -jav indicates the average over the quenched 
disorder. 

Two different routes were used to check numerically the above scenario in three 
dimensions: One way is to perform Monte Carlo simulations of finite but large 
systems above Tc, where the correlation length ^ is smaller than the system size 
L, and to perform a scaling analysis of the spatial correlations G{r). This has 
been done by Ogielski and Morgenstern | J208| | in extensive simulations on a special 



purpose computer. The other way is to confine oneself to smaller system sizes and to 
apply finite size scaling methods to analyze the data for the spin glass susceptibility 
XsG and certain combination of its moments, which has been done by Bhatt and 



Young |2^, Q . Both approaches come to the conclusion that there is indeed a finite 
temperature spin glass transition in the three-dimensional Ising spin glass model (|I]) 
with a binary distribution Jij = ±1, the critical temperature is Tc = 1.175 ± 0.025, 
the correlation length exponent is z/ = 1.3 ±0.1 and the exponent for the spin glass 
susceptibility is 7 = 2.9 ± 0.3. The results from high-temperature series expansions 



||270| , |271| , P94] , |153|| agree with these values within the errorbars. 



At a second order phase transition the characteristic time scale r for the decay 
of on-site correlations 

qit) = [{Si{t)S^iO))] ~ r^qit/r) with x = (d - 2 + 7])/2z (4) 

is also expected to diverge according to r ~ ^^ ~ (T — Tc)~^'^ with z being the 
dynamical exponent. One possibility of defining an effective relaxation time is r = 



U^dttq{t))/{!^dtq{t)), for which Ogielski [gOg] finds z = 6.1 ± 0.3 for the 3d 
±J spin glass model. This is an unusually large value for z (which causes serious 
equilibration problems in Monte Carlo studies of the critical properties of this model) 
and an exponential divergence of r at Tc like r ~ exp{A/(T — T^} cannot be 
excluded. Bernardi and Campbell [^, [l^ perfomed a similar analysis for other bond 



distributions and found differences for the value of the exponent x depending on the 
particular form of these distributions, indicating that universality with respect to 
the latter possibly does not hold. Hukushima and Nemoto ||127|| investigated with 



the same method the two-dimensional case, where an exponential divergence of r at 
T = was found. 

In a more recent work Blundell et al. |^6|] perform a scaling analysis of the non- 



equilibrium critical dynamics based on the hypothesis that the equal time spin-spin 



correlation function measured at time t after a temperature quench from T = oo to 
T = Tc obeys 

C{r,t) = [{S,it)S.^rit))X. = r-^'-'^^^fir/t^'^) . (5) 

Note that here (■ ■ ■) means an expectation value with respect to the time dependent 
probability distribution generated by the underlying stochastic process and not the 
thermodynamic equilibrium expectation value. With other words, according to the 
hypothesis (^ one does not need to equilibrate the system to get estimates for 
the equilibrium critical exponents t] and z. In this way a value z = 5.85 ± 0.3 and 
T] = — 0.29±0.07 was obtained, concurring with the value reported above by Ogielski 
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Another approach, suggested by Bhatt and Young [|^ , is to look at certain ratios 



of moments of time-dependent spin autocorrelation functions as for instance 



av 



Since this is a dimensionless quantity it is expected to scale like R'{t) = R{t/T) with a 
characteristic time scale r = t{L,T) oc L^f{L^/'^{T — Tc)). Again the determination 
of equilibrium critical exponents proceeds via these hypothetical scaling forms of 



non-equilibrium quantities. Assuming that T^. = 1.175 (see above) in 1 25 1 it is found 
that z = 6.0 ± 0.5, again in good agreement with Ogielski's result [p09|| . 

As we have seen so far there seems to be rather compelling evidence that the 
three-dimensional Ising spin glass in vanishing external field has a second order 
phase transition at a finite temperature T^. However, a recent Monte Carlo study 
by Marinari et al. [|172|| of the three-dimensional Ising spin glass on a body-centered 
cubic lattice, where second and third nearest neighbor interactions have been taken 
into account in order two increase the hypothetical critical temperature, revealed a 
less clear picture: It was found that in this model the data are compatible with a non- 
zero temperature phase transition (with T^ = 3.25 ± 0.05, which is a non-universal 
quantity and therefore different from the value reported above, z/ = 1.20 ± 0.04 and 
7 = 2.43 ± 0.05) as well as with a zero temperature transition: 



XSG - A[exp{B/TY-l]+C 
^ ~ a[exp(6/T)^-l]+c. 
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Marinari et al. | 172 | report that this form fits their data better with a smaller 



systematic error, thus giving new support for the speculation that three dimensions 



may well be (instead of close to but larger than) the lower critical dimension of the 
short ranged Ising spin glass model. 

In any other number of dimensions the situation is much clearer: In two dimen- 
sions it was shown that no phase transition is present at any finite temperature 
Tn\ , ITTDj §§ gUDllTSg, ||, [T^, |2S3|, |2g. For a ±1 distribution the most recent re- 
sults are from Bhatt and Young \H\ reporting that Tc = 0, u = 2.6 ± 0.4 (^ ~ T^'^) 



and 7 = 4.6 ±0.5 {xsg ~ T^"/). A continuous distribution yields significantly differ- 
ent results indicating that this case is in a different universality class than the ±1 
distribution, since the latter has a large ground state degeneracy. In four dimensions 
one finds much stronger indications for the existence of a phase transition than in 

For a Gaussian distribution it is found that Tc = 1.75 ± 0.05, 
1.8 ± 0.4 ^ and z = A.8 ± 0.4 [|^. Furthermore the ex- 



three dimensions [^ 

u = 0.8 ±0.15, 7 = 

tensive investigation of the probability distribution P{q) of the replica overlap q 

||226| , |220| , |63[| revealed a picture that is reminiscent of the Sherrington-Kirkpatrick 

(SK) model [|262|| and its solution with broken replica symmetry ||181|| . The upper 

critical dimension of the Ising spin glass, above which mean field results should ap- 



ply, is expected to be six [|117]| . This prediction has been confirmed by Monte Carlo 
simulations of the six-dimensional Ising spin glass with ±1 couplings ||295|| , where at 
Tc = 3.03 ± 0.01 mean field critical behavior is found, which is i/ = 1/2 and 7 = 1, 
with logarithmic corrections due to being at the upper critical dimension. 

The effect of an external magnetic field on the critical behavior of the Ising spin 
glass is also far from clarified. The mean field scenario [ p.81|| would suggest the 
existence of an AT (de Almeida- Thouless |^) line in the /i-T-phase diagram and 
thus a spin glass transition even within a non- vanishing external field. On the other 
hand phenomenological models ||89|, ^] dispute the relevance of mean-field results 



for finite dimensional spin glasses and predict the absence of a phase transition in 
a field. All numerical results obtained so far for four dimensions [p.l0| , P, ^ |240| 



seem to be compatible with the existence of a spin glass transition at a nonvanishing 
temperature in a low enough external field. 



In three dimensions the situation is still unclear: Kawashima et al. ||144| , p.45|| have 
investigated the variance of the probability distribution of the central spin magne- 
tization as well as the probability distribution of the replica overlap and observe 
that for increasing system size it scales to zero in non-vanishing external fields. This 
implies the non-existence of a phase transition. Following a proposition by Singh 
and Huse ||272|| , Grannan and Hetzel ||110|| looked at the contour lines Xsg{T, h) = 1 
in the T-h diagram and observe an upwards bend (i.e. h -^ oo) for T — i> in four 



dimensions (which indicates a spin glass phase) and a downwards bend (i.e. h ^ 0) 
in two dimensions (meaning the absence of a spin glass phase). This is depicted in 
figure |l|. The three-dimensional case is only marginally bending upward, suggesting 
that three is close to but larger than the lower critical dimension also in an external 
field (see also ||120|| for the same study of the 3d Ising spin glass on the face-centered 
cubic lattice). Ritort [25C] came to simlar conclusions by investigating the quanti- 



ties describing the sensitivity of the spin glass phase with respect to magnetic-field 
perturbations. Finally Caracciolo et al. ||^, ^ Q presented results of Monte Carlo 
simulations of the (three-dimensional) system 



which is identical to the two-replica method used by Bhatt and Young ||22| apart 
from the existence of a coupling of strength e between the two replicated systems a 
and T. This parameter e plays a similar role as the magnetic field in a ferromagnet 
— it is the conjugate field to the order parameter q, the replica overlap. Hence 
one expects a discontinuity in q{e) in the limit e ^ for temperatures T < T^, if 
there is a phase transition. The authors conclude that their results are inconsistent 
with the predictions of the droplet model [^ and support the mean field picture 
including replica symmetry breaking and the existence of an AT-line (i.e. a spin glass 
transition within a field) even in three dimensions. These claims have been heavily 
disputed by Huse and Fisher ||132|| (see also the subsequent reply by Caracciolo et 
al. ^1), the main counter-argument being that their simulations are performed at 



temperatures above the hypothetical AT-line and that the system sizes are much too 
small to give reliable results for the low field regime. We will return to the issue of 
the correctness of the droplet theory in finite dimensional spin glasses below. 

3 Ising spin glasses: Non-equilibrium dynamics 

Most experiments on spin glasses at low temperatures are performed in a non- 
equilibrium situation due to the astronomically large equilibration times. Since the 
seminal work of Lundgren et al. [|169|| it has been realized by the experimental- 



ists [|170| , |126|| that magnetic properties of spin glasses strongly depend on the time 
they spent below the spin glass transition temperature Tc, a phenomenon that has 
been dubbed aging, (see e.g. [|164| , p92|| for a review about the experimental situ- 



ation). Figure |^ shows the result of a typical aging experiment: The spin glass is 
rapidly cooled to low temperatures in an external field h (in three dimensional spin 
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glasses below the transition temperature) and the field is switched off only after a 
macroscopic (several hours) waiting time t^. As soon as the field is zero (which de- 
fines the time t = 0) the time-(t)-dependence of the thermoremanent magnetization 
-/^trm(^,^«)) is measured. Usually (as for instance in simple ferromagnets) the equi- 
libration time is microscopic and for macroscopic waiting times one would measure 
always the same magnetization curve Mtrm(^, iw) independent of t^. As can be seen 
in figure ^ this is completely different in spin glasses, which is an obvious manifes- 
tation of the huge time scales of the glassy dynamics. In fact a spin glass transition 
is not a necessary ingredient for this scenario to occur, it has also been reported 
in experiments with two-dimensional spin glasses like Ru2Cuo.89Coo.11F4 |^ ^ or 
CuMn- films [|175|| and in models without frustration [|238|| or disorder ||15CI| , |73| . The 
observation characteristic for aging can also be made in other amorphous or strongly 
disordered materials like polymers [p8CI|| , charge-density- wave systems []26| , |2^ or 
high-temperature superconductors ||241 . 



In Monte Carlo simulations the spin autocorrelation function 



C(t, ty 



1 



N 



Ni 



Y.[{s,{t + t^)s,{t^))l 



(9) 



i=l 



is the quantity that is most convenient to investigate the non-equilibrium dynamics 
of spin systems. Here, as in equation (|), (■ ■ ■) formally means an average over 
various thermal histories and initial conditions. However, as long as the number of 
samples for the disorder average [■ ■ -jav is large enough, it does not make a difference 
(numerically) if one neglects this thermal average. 

The limit t^ — * 0, which corresponds to experimental measurements of the re- 
manent magnetization Mtrm(''^) after removing a strong external field (of the order 
of the saturation field strength), measures the configurational overlap with a fully 



magnetized initial state. Early investigators of this quantity ||2^, |146| , p.47|| observed 
already an algebraic decay in a certain time window at low temperatures. This has 
been confirmed recently for the SK-model [B^, 218 | and finite dimensional spin glass 
5^, 218 |. In more extensive simulations for the 3d EA-model | 234 | the 



models 



temperature dependent exponent A(T) in 

MTRM(t) ~ t-^(^) , (10) 

has been determined for various temperatures below the spin glass transition temper- 
ature and found a good agreement with those reported for the amorphous metallic 
spin glass (Fea;Ni(i„^))75Pi6B6Al3 [ p.07|| , which is particularly well suited for measure- 



ments of the remanent magnetization also for low temperature since its saturation 
field strength is reasonably low (see figure ^. 
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For non-vanishing waiting times t^ 7^ the function (^ shows a behavior that 
is characteristic for aging phenomena. Andersson et al. and Rieger p34|| showed 
that quahtative features of experimental observations can be reproduced by Monte 
Carlo simulations of finite dimensional spin glass models. In addition to C{t, 1^) the 
magnetic response function x(t,t^) = Mtrm(^,^w)/^ has been determined, where 
A^trm(^, tw) is defined in the introductory paragraph of this section. Its logarithmic 
derivative 9x(t, t^)/(91og t possesses a maximum at t = t^, as observed in the cor- 
responding experiments ||199| , |1(J8| | . Furthermore one observes that the fluctuation- 
dissipation theorem xi^jtw) = P{^ — C{t,tw)) is violated for t > t^. A more quan- 
titative analysis of the the data for the function C{t,tw), of which a typical set is 
depicted in figure §, reveals the following picture ||234|| : In three dimensions there 
are strong indications for the scaling relation 

C(t,tJ=t-^(^)$T(t/tJ with $t(2/) ~ I ^X'aV) III IZl^ (11) 

to hold This characteristic t/t^^j-scaling has also been observed in experiments with 



the insulating spin glass CdCr1.7Ino.3S4 [|206| , g], in two-dimensional spin glasses ||237 



the SK- model |7l|], in simplified spin glass models pO| , |219|| and one-dimensional 
models |p38|| . For a spin glass below the critical temperature the relation ([TI|) is 
expected to hold for any finite waiting time t^ (as long as the system is infinite), 
which expresses the fact that for all temperatures below Tc the equilibration time is 
infinite — in contrast to the situation in e.g. simple ferromagnets. For t^ ^ 00 one 
obtains the equilibrium autocorrelation function limi^_^oo C'(t, t^,) = q{t) ~ t~^^'^\ 
see equation (^). It should be noted that in spin glasses with replica symmetry 
breaking a more complicated scenario can occur, see ||9^ for a discussion. At Tc one 
has x{Tc) ~ 0.07, which equals the exponent reported for the equilibrium dynamics 



of the 3d EA-model ||209|| as well as from experiments on the short-range Ising spin 



glass Feo.sMuo.sTiOs [ p.l4|| . In other words: as long as t -C t^ one observes a quasi- 
equilibrium dynamics (on time scales smaller than the waiting time t^) and the 
crossover at t = t^„ signals the onset of the non-equilibrium behavior characterized 
by an exponent A(T) that is significantly larger than the corresponding equilibrium 
exponent x(Tc) (e.g. A(Tc) ~ 0.39, for this value see also [p.29|] ). 



At this point a short excursion on semantics might be appropriate. A: When 
speaking of equilibration time in the preceeding paragraph we mean the equilibration 
time within one ergodic component of an infinite system, like the one with positive 
(or negative) expectation value for the magnetization in a ferromagnet below Tc. 
This definition is independent of the answer to the question, whether the phase 



space of a finite dimensional spin glass splits into an infinite number of ergodic 
components separated by infinite free-energy barriers below T^ (as the SK-model) 
or not. We observe that \im.t-*ooC(t,tw) -^ for any waiting time t^, but the 
more natural expectation for dynamics within an ergodic component would be an 
exponential decay of q{t) to a finite value qea- However, no indication for this has 
been reported yet [gU^, |2^, |7T1. 

B: The scenario ( ]TT[ ) is what we call aging, and we use the expression interrupted 
aging for a situation in which we have some finite equilibration time Tgq in such a 
way that ([iTl) holds approximately as long as t^ -C Teq and is simply replaced by 
true equilibrium dynamics C{t,t^) ^ q(t) for t^ ^ r^q. This is what happens in 
two-dimesnional spin glasses ||247| , |248| , p37|] or random ferromagnets ||238|| , but also 



in pure systems like the one considered in ||150|| , where depending on the system 



parameters the time r^q can be extremely large. Obviously for an astronomically 
large Teq neither experiments nor Monte Carlo simulations might be able to discrim- 



inate between aging and interrupted aging. It has also been noted |[L5CI|| that in the 
pure ferromagnetic Ising chain the aging scenario ([TT|) holds at T = with a; = 0, 
thus aging phenomena do not rely upon randomness or frustration althought the 
latter can greatly enhance it. We would like to stick to this nomenclature in this 
review. However, there exists also a different perception of what aging might be 



|[70| , [7^ , p.83| , p6[] , see also the discussion following eq. ( ]T6D below. In essence, aging 



understood in this way is a property of off-equilibrium dynamics that becomes man- 
ifest only in particularly chosen limiting procedure for the times t and ty^. A scenario 
like ([TTl) with x(T) > would be called "interrupted aging", even if Teq = oo. This 
concept of aging is motivated by analytically tractable mean-field models, where 
various infinite time limits can be done in a straightforward manner — its relevance 
for real spin glasses and finite time experiments remains to be tested. 

The underlying mechanism for the aging scenario ([TI|) in finite dimensional spin 
glasses is a slow domain growth, as suggested by Fisher and Huse [^ proposing (ad 
hoc) a logarithmic growth law for the characteristic domain size R{t) oc (logt)^/"^. 



and by Koper and Hilhorst ||157|| , stipulating (ad hoc) an algebraic growth R{t) oc 
f^C^), Preliminary investigations on this issue [|131|| find some numerical support 
for the first hypothesis, however, the above mentioned facts (the asymptotically 
algebraic decay of C(t, t^) and the t/tw scaling) and recent very extensive simulations 
pledge more in favor of an algebraic growth: Rieger et al. [ p37|| perform Monte 



Carlo simulations of the two-dimensional EA-model with Gaussian couplings and 



investigate the growth of spatial correlations via the correlation function 

1 ^ 
G{r,t^) = J;,Y.l(S^{tu,)S^+rK))X. , (12) 

which becomes G{r) of equation (0) in the limit t^ — ^ oo. Note that the domain 
growth in various strongly disordered systems like the site-diluted Ising model, the 
random field Ising model and the random bond ferromagnetic Ising model, has been 
investigated frequently via Monte Carlo simulations (see the review by Chowdhury 
and Biswal in this series [§31 )• These models have the advantage that their ground 
state is known to be ferromagnetic, which makes the identification of domains easy. 
However, this is not the case for the present system and such an investigation is much 
more difficult here. Hence, a straightforward way to determine a typical length scale 
for spatial correlations in a spin glass is to calculate the function G{r, t^) defined in 

(0)- 

In order to take into account the square of the thermal average in (|12]) in [p37|| two 

replicas a and b of the system were simulated. Then [(S*," (t^)S'f (t^)5'f_,_r(ti„)5'f^r(^w))]ai 
instead of [(5'j(tu;)5'j+r)^]av was calculated, giving the same results with a much bet- 
ter statistics for the first quantity. To improve the statistics one also has to average 
over a suitably chosen time window around t^. The correlation length is defined via 
^(tw) = 2 /q°° drG{r, tw). It turns out (see figure ^ that 

e(t.) oc r(^) , (13) 

with an exponent a{T) that decreases with temperature yields a good fit to the 
time dependency of the correlation length. Of course in the two-dimensional spin 
glass, which does not have a phase transition at any finite temperature (see previous 
section), ^(t^) has to saturate at a finite value. This, however, happens at small 
temperatures (T ~ 0.2) only at very large times t^ inaccessible to computer simu- 
lations. It should be noted that a logarithmic fit ^{tw) ~ (logt)^/'^ also works fairly 
well — with a temperature independent exponent ip ~ 0.63. In recent experiments 
on the two-dimensional spin glass Ru2Cuo.89Coo.11F4 |p47| , p4^] the scaling behavior 



of the ac-susceptibility has been analyzed and it turned out that this is compatible 
with a prediction made by Fisher and Huse |^ assuming a logarithmic growth with 



ip ^ 1.0. Unfortunately no direct measurements of the correlation length itself are 
available experimentally up to now. 

For three dimensions the same kind of calculations as described above have been 
done lp.5111 (see also the work by Sibani and Andersson ||268| , || discussed below) 



which indicate the validity of algebraic growth law (|T3D, too. A continuous set of 
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exponents like <y{T) for domain growth (or sometimes called coarsening) might be 
disturbing especially if one considers the work of Lai, Mazenko and Ma ||163|| . They 
classify various non-equilibrium systems according to the underlying scaling law for 
(free) energy barriers B{L) that domains of size L{t) have to overcome in order to 
grow further: two of them lead to L{t) oc ^/t and two lead to logarithmic growth — 
the fourth being identical to the proposal of Fisher and Huse [0] for spin glasses 
B{L) ~ L'^, which leads to L{t) ~ (logtY^'^ via activated dynamics. However, as 
suggested by Rieger ||234|| , the scaling law 



B{L) = A(T)logL 



(14) 



leads (via an activated dynamics scenario in which it takes a time r ~ exp{B/T) 
to overcome a (free) energy barrier B) to an algebraic domain growth with a con- 
tinuous set of exponents as in (|l^). Moreover, following [^ in their argumenta- 
tion, we stipulate for the remanent magnetization that Mtrm(^) ~ L{t)~^ and for 
the autocorrelation-function that C{t,tu,) ~ [L(t^t,)/L(t + t^)]^. Then it follows 
Mtrm('^) ~ t^'^^'^\ C{t,tw) ~ t~'^(^) for t ^ t^ and also the t/tw scaling stated in 
(p]l]). The following table gives a collection of predictions for the different scaling 
assumptions. 





Droplet-model |91 


MC-sim. 


Energy barrier 


B 


AL^ 


A(T)logL 


• 


Activated dynamics 


r 


T ~ exp B/T 


r ~ exp B/T 




Domain size 


R(t) 


(fi°g')"* 


^a(T) 


(•) 


Remanent Magnetization 


^TRM 


(f>°.0"'* 


t-HT) 




Aging 


C(t,t^) 


Ty / log i+i„\ 
^ V logt™ J 


^(i) 




Asymptotic decay 


t>t^ 


(log t)-'/^ 


^-A(T) 






t<t^ 


(logt)-^/^ 


^~x{T) 


(•) 



The bullets indicate that the Monte Carlo simulations described above and below 
confirm the corresponding prediction in the last colomn. A bullet with brackets 
means that the numerical data can also be interpreted according to the predictions 
of the droplet model. Concluding the hypothesis (|I^) seems to give a more consistent 
description of experimental and numerical results on aging in spin glasses presented 
so far. 

Recently Sibani and Andersson [^6^ , ^ gave further support for the relation 
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(0) by means of the following procedure: In Monte Carlo simulations of the two- 
and three-dimensional spin glass reference states \E' are generated with the help of a 
careful annealing procedure down to T = 0. Then a certain amount of spins is flipped 
randomly under the constraint that the energy difference of the new state to the old 
reference state remains below a certain lid b. Via zero-temperature dynamics the 
system has the opportunity to relax into a local energy minimum configuration that 
is stable against single spin flips. Finally the new state is analyzed by identifying 
all clusters of reversed spins. For these clusters the size and energy distribution 
is recorded. It turns out that the cluster volume increases exponentially (or even 
superexponentially) with the lid b, implying a logarithmic dependence for the energy 
barriers as in ([T^ . This dependency disagrees with the result of Gawron et al. 



lOOJI , who use an exact search algorithm for calculating the energy barriers against 
inversion of ground states. They obtain in two dimensions B{L) ~ L, which means 
ip = 1. A serious caveat in their approach might be the fact that their investigations 
have been constrained to system sizes L < 5. 

Especially with regards to the latter investigations it seems that the system- 
atic analysis of the phase-space structure of short range spin glasses, especially their 
ground states and low lying excitation, seems to become feasible with increasing com- 



puter power These studies already gave valuable insights ||266| , |267| , |154| , |155|| , sup- 
porting the phenomenological theories of hierarchical relaxation for the spin glass dy- 



namics a la Schreckenberg p49|| , Ogielski and Stein |p07|| , which have been improved 



further by Sibani and Hoffmann [ p65| , p51| , |124|| in order to model simple aging, tem- 
perature step experiments and violation of the fluctuation dissipation theorem. In 
particular with regards to improved optimization algorithms [|1 1 1| , |lO| , |281| , |269|] (see 
also ||143|] and the articles by Stariolo and Tsallis [p76|] , Hogg ||125|] and Tomassini 



288|] in this series) significant progress in the investigation of ground state properties 



in spin glasses can be expected in the future. 

Another procedure devised to test various phenomenological spin glass theories 
are so-called temperature cycling experiments. They consist of two temperature 



changes during the time in which the material is aged in the spin glass phase |225 



either a short heat pulse is applied to the spin glass during the waiting time after 
which the relaxation of e.g. the thermo-remanent magnetization is measured, or a 
short negative temperature cycle is performed, which is the same as a heat pulse 
but with a negative temperature shift during the pulse. It has been pointed out 
| 166|| that this kind of experiments can discriminate between the droplet picture |91 



and the hierarchical picture [p.64|| (see also ||134|| and [Q for a microscopic theorie 
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of adiabatic cooling, which are also capable to explain some of the obsevations) . 
The interpretation of the experimental situation is still controversial (see also ||296|| 
for a refreshing description of the present situation) — essentially the Saclay group 
pledging in favor of a hierarchical interpretation ||225 , 166 , p92|| and the Uppsala 



group taking the part of the droplet model ||109| , |175|| . The situation is the same 



with regards to numerical simulations of temperature cycling experiments: Rieger 
1 235] comes to similar conclusions as the former and Andersson et al. ^ interpret 



their results in full agreement with the latter experiments. In both Monte Carlo 
simulations the same quantities are studied (for instance the thermoremanent mag- 
netization Mtrm(^; tage), whcrc tage HOW meaus the whole time in which the temper- 
ature cycle has been performed), but the resulting data are analyzed differently. An 
alternative approach to the study of the temperature dependence of the spin glass 



state has been suggested by Jan and Ray ||140|| using damage spreading (for the use 



of the latter concept in Ising spin glasses see [^, 0, |5^ and the article by Jan and 



de Arcangelis in this series [|141|| ). 

It might be difficult to make conclusive progress in this direction, because these 
kind of experiments (real and numerical) always yield ambiguous results that invite 
to one or the other interpretation. However, the predictions of the droplet theory 
heavily rely on the concept of chaos in spin glasses ^\, [L97| , |198|| (meaning the 



significant sensibility of the spin glass state to infinitesimal changes of parameters 
like temperature or field) that can be quantified by an overlap length ({T,AT) via 
the hypothesis 

[{SiSi+r)T{SiSi+r)T±AT]av ~ exp{-r/C(T, ±AT)} , (15) 

where (■ ■ ■)t means the equilibrium expectation value of one system at temperature 
T. The characteristic length scale C(T, ±AT) should be finite for the positive and 
negative sign of AT and it should decrease with increasing AT. There are quan- 
titative predictions for these dependencies, inaccessible to experiments, which can, 
however, be calculated in Monte Carlo simulations. This seems to be a promising 
endeavor. It should be noted that the existence of an overlap-length as defined 
in ([T5|) is a prediction of phenomenological concepts [H^, RH, |91I] and microscopic 



theories on the basis of the SK-model, as has been reported by Kondor ||156|| and 



discussed further by Ritort |240 



It should be mentioned that very recently very appealing theoretical concepts 
for the non-equilibrium dynamics in spin glasses have been developed. The phe- 
nomenological theory that fits the experimental data for aging experiments in the 
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best way up to now was introduced by Bouchaud [HOl O, E3[. In essence it is 



a 



diffusion model in an abstract space in which each state is characterized by a ran- 
dom (free) energy and hence by a random, exponentially distributed trapping time. 
By arranging them in a tree like structure, very reminiscent to the phase space 



structure that emerges from Parisi's solution of the SK- model [TEl], one obtains 
functional forms for MTRAtit, tw) and C(t, t^) that yield excellent fits to experimen- 
tally measured data. Schreckenberg and Rieger ||25(J| propose a different diffusion 



model, which is based on an ultrametric tree that incorporates a separation between 
quasi-equilibrium and non-equilibrium branches. In this way aging occurs natu- 
rally via exploration of quasi-equilibrium sub-branches of increasing depth. Also 
the microscopic theory for the off-equilibrium dynamics of mean field models of spin 
glasses has been pushed forward by Cugliandolo and Kurchan [^, ^ and Franz and 
Mezard ||183| , [96[| : The mathematical difficulties for an analytically exact solution of 



the dynamical off-equilibrium mean field equations for e.g. the SK-model come from 
the lack of the fluctuation-dissipation theorem (FDT) that relates autocorrelation 
and response function (which allows the analytical solution in case of equilibrium 
dynamics ||274| , |273| , [68| , |148|| ). The new approaches circumvent this by considering 



a so-called fluctuation-dissipation ratio defined via 

and postulating a particular set of properties for this function x{t,t') in various 
asymptotic limits, essentially setting up an "ultrametric" for timescales. In this 
way Parisi's static, equilibrium (!) order parameter function q{x) finds its counter- 
part in off-equilibrium dynamics. It will be a challenging endeavor to test these 
interesting ideas via Monte Carlo simulations and to check, whether they might also 
be applicable in finite dimensions. The analysis of Monte Carlo results for the three- 
dimensional EA spin glass are compatible with the proposed scenario [^ although. 



due to the marginality of the three-dimensional case, further investigation in higher 
dimensions are desirable and would certainly give a much clearer picture. 

Finally we would like to point out that aging and glassy dynamics have gained 
much interest very recently. Apart from the issues mentioned above in connec- 
tion with spin glasses the central point of the research activities is to model glassy 
behavior with spin systems that have no quenched disorder (in order to try to un- 
derstand the glass transition that leads for instance to the formation of window 
glass): among them are two- and three-dimensional models with competing nearest 
and next-nearest neighbor interactions ||259| , |264|| , the anisotropic kagome antifer- 
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ro magnet | 229 , |58| , one-dimensional spin models with p-spin interactions | 150 |, the 
Bernasoni model (also being a p-spin interaction model, but with infinite range in- 
teractions) [^, |173| , |185| , 159 1 and certain field theoretial models with infinite range 



interactions [73, 98]. Apart from the latter two, the bulk of investigations of this 



very fascinating subject has to be done numerically via Monte Carlo simulations. 

4 Numerical recipes 

It is obvious that thorough Monte Carlo simulations of spin glasses require a huge 
amount of computer time. As long as one is interested in purely dynamical phenom- 
ena as described in the previous section, there is no other way than to implement 
some algorithm that generates a stochastic process described by the Master equation 
for the probability distribution of the spins. The possibility that is most frequently 
used to achieve this, is the heat bath, Glauber or Metropolis algorithm. Here a 
random number generator is used to decide, whether a certain spin configuration S_ 
is modified or not, according to predefined transition probabilities, as for instance 

w{S -^ S!) = min {1, exp(-/3AE)} , (17) 

where AE is the energy difference between the old {S_) and new (S') configuration. 
For Ising systems very efficient implementations of this algorithm with single spin 
fiips have been presented that reach a speed of 3 ■ 10^ spin update attempts per 
second on a single processor of a Cray YMP: for the pure Ising model by Ito and 
Kanada [ |13711 and Heuer [ |121|1 , for the random field Ising model [|231|| and finally for 



a whole class of Ising models with or without quenched disorder defined via binary 
variables ||138|| (see also the book of de Olivera ||213|| for a review on multi-spin 



coding techniques). The enormous speed up compared to older multi-spin coding 
techniques is based upon the use of a single random number for different (up to 64) 
systems, like e.g. 64 different samples of random bond configurations of the EA spin 
glass. This method does not work for a continuous probability distribution of the 
quenched disorder variable (like Gaussian) — in this case more conventional methods 
have to be used, which are at least one order of magnitude slower. Although nearly 
all results reviewed above have been obtained with these more or less optimized 
algorithms they reveal serious deficiencies at low temperatures: 

Equation (|1^) means that the new configuration S! is rejected if a random num- 
ber, equally distributed between and 1, is larger than w{S_—^ 5'), therefore this is 
called a rejection method. In system with a huge number of local energy minimum 



15 



configurations the number of rejections becomes very large at low temperatures and 
most of the CPU time is wasted by generating random numbers and calculating (or 
looking up in a table) the transition rates. The only effect is to increment the time 
that has passed since the beginning of the simulation by one unit. Bortz, Kalos 
and Lebowitz |^, ^ suggested already 20 years ago a method that seems to be 
more efficient: They proposed to accept a new configuration always and then to 
increment the time by a random number At obeying a probability distribution that 
is characterized by the sum over all probabilities for transitions away from the old 
configuration: 

P(At) = T-^ exp(-At/r) with t'^ = ^ w(5' -^ 5') (18) 

s! 

Although this idea might be capable of circumventing the above mentioned slowing 
down it has not been used very frequently up to now [KI^, 264, ^ 16(J| ]. One of 



the problems that typically occur also here is that once a transition away from a 
local minimum configuration has taken place the system will very soon return to it 
and thus decrementing the efficiency of the algorithm. Very recently Krauth and 
Pulchery [|159|| proposed a variant of this method that seems even to avoid this 



caveat of short cycles: By keeping track of the configurations already visited (which 
means storing them into the computer memory) the algorithm is forced to generate 
new configurations in each iteration. The memory is cleared as soon as a lower 
local energy minimum is encountered. In this way they were able to explore times 
scales equivalent to 10^^ conventional Monte Carlo steps for a particular spin model 
with 400 spins. This is very promising and would mean a real breakthrough if such 
a performance could also be achieved with finite dimensional spin glass models of 
reasonable size. 

A completely different method that tries to avoid the tremendous slowing down 
caused by coexisting energy minima in phase space with large energy barriers be- 
tween them is the so-called multicanonical ensemble |T^ or simulated tempering 



| 171|] , which is suited for the calculation of equilibrium properties and also ground 



states of spin glasses. Usually (as in the above mentioned algorithms) one performs 
an importance sampling with the canonical distribution Pcan.(-E') ~ exp{—E/T), 
which is sharply peaked around one (at high temperatures T), two (in case of first 
order phase transitions) or several (in spin glasses or random field models) values 
of the energy E. In order to escape from one energy minimum region to explore 
others and, in particular for interface problems, the regions in between them, one 
tries to generate an ensemble Pm{E) that is approximately fiat for the energy inter- 

16 



val of interest, especially between the sharp maxima of the canonical distribution. 
This can be achieved with appropriately chosen functions a{E) and (3{E) in the 
multicanonical ensemble 



P^{E)r^exp[-f3{E)E + a{E)] 



(19) 



Since they are unknown a priori, they have to be determined recursively during the 
simulations \TE, TD]. Finally, by re-weighting with exp[—E/T + (3{E)E — «(-E)], 
the canonical distribution can be reconstructed. With regards to spin glasses, this 
ensemble has proven to be useful in particular for the investigation of ground state 
properties [ll9| , ^ |2l|| (for this problem see also the above mentioned works ||111 
ll. Hi], m, |27|, ini, Hi])- 



Cluster algorithms as suggested by Swendsen and Wang 284|| have become very 



useful in the investigation of pure systems (see [p85|| for an overview). However, they 
do not yet give a significant improvement over local algorithms simulating single spin 
flip dynamics in spin glass or random fleld models — with at least one exception: 
In a special implementation of a cluster algorithm for the two-dimensional Gaussian 
EA-model Liang showed ||167|| that the logarithm of the relaxation time is flve times 
smaller than the usual Metropolis algorithm. But he also pointed out that this 
efficiency will not be reached in higher dimensions. A cluster algorithm for the 
random fleld Ising model has also been proposed ||81|, but its efficiency remains to 
be tested. The reason for this is that due to randomness and frustration it is not 
obvious at all how to construct spin clusters that can be reversed with acceptable 
rates. Since the correct construction of these clusters is the main problem in spin 
glasses we should also mention the work on various cluster concepts in disordered 
systems [|7^, |2|, |27|, [m ]. 

Finally let us mention histogram techniques ISH] that also have been used fre- 



quently for pure systems |p85|| : in principle they allow to get information about 
thermodynamic quantities of one system (i.e. one realization of the quenched disor- 
der) in a whole temperature interval by a Monte Carlo run at one single temperature. 
However, this needs a lot of book keeping already for pure systems, where several 
thousand configurations have to be stored. In spin glasses and random field systems 
in addition one has to perform an average also over this disorder, which might cause 
a serious difficulty for data storage. 
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5 Random field systems 

If one mixes a typical antiferromagnet like FeCl2 with an isostructural nonmagnetic 
material like C0CI2 or NiCl2 one obtains diluted antiferromagnets, which are dis- 
cussed in this series p2|, [^57||. Within an uniform external field, however, a large 



degree of frustration is induced and a completely new universality class (with re- 
gards to critical properties) emerges. Fishman and Aharony |9^ and also Cardy 



||55|| pointed out that it is the same as the universality class of the random field Ising 
model (RFIM), which is defined by the Hamiltonian 

H = -JY.SiSj-Y.h,S,, (20) 

where the first sum is again over nearest neighbor pairs on a d-dimensional lattice 
and the hi are independent quenched random variables with [/ii]av = and [hf]av = 



h^. A simple heuristic argument by Imry and Ma [ p.33|| shows that the lower critical 
dimension of the RFIM (^) should be li = 2. Indeed it has been proved rigorously 
^ that in three dimensions for small enough random field strength h there is 
ferromagnetic long range order at low enough temperatures. Thus the existence of a 
phase transition in three (and higher) dimensions is assured (in so far the situation 
is slightly better than in spin glasses), however, there has been a long lasting debate 
on the critical properties, which still awaits a solution. 



Villain ||291|| and Fisher [^ proposed a scaling theory for the random field tran- 
sition that relies upon the assumption that random field induced fluctuations dom- 
inate over thermal fluctuations at T^. This implies for the singular part of the free 
energy 

^sing ~ e' , (21) 

where ^ is the correlation length and 6' is a new exponent. Random field fluctuations 
alone produce typically an excess field of the order of ^'^/^ within a correlation 
volume, so a naive guess would be 9 roughly 1.5 in three dimensions. From (pT]) the 
modified hyperscaling relation follows 

2-a = v{d-e), (22) 

and it also implies an exponential divergence of the relaxation time r at Tc- t ~ 
exp(A/^''). The decay of the connected ([(5o'^r)]av - [(^o)(^r)]av ~ r-('^-2+'')) and 
disconnected ([(>S'o)(S'j.)]av ~ r"^'^"^"'"'')) correlation functions at Tc define two expo- 
nents ?7 and rj. These are expected [p91j , |88[| to be related to the new exponent 9 



via 

e = 2-r] + ri . (23) 



And, as usual, one has the scahng relations 7 = u{2 — rj), P = (rf — 4 + fj) and 
a + 2/5 + 7 = 2. Obviously there seem to be three independent critical exponents 
and a central issue of the activities on the critical properties of the RFIM is the 
quest for an additional scaling relation. Already Imry and Ma ||133|| conjectured 
that Fsing ~ X? where x ~ (T — Tc)^'^ is the susceptibility. This would imply 
e = 2-r] {^ and therefore via (^) 

fi = 2r] . (24) 



A set of analytical arguments by Schwartz et al. |p52| , |253| , |254| , |255|| supports this 
two-exponent scaling scenario indicated by p^). Hence the exact Schwartz-Soffer 
inequality ||253|| fj < 2ri might be fulfilled as an equality. Results of numerical 



studies are compatible with this picture: In Monte Carlo simulations of the the 
diluted antiferromagnet in a uniform field (DAFF) Ogielski and Huse [plO|] found 
1] = 0.5±0.1 and Tj = 1.0±0.3. Ogielski [|211|| investigated ground state properties of 



the RFIM via combinatorial optimization methods |^ and found rj = 1.1 ± 0.1 and 
z/ = 1.0 ± 0.1. Rieger and Young [ p32|] performed the most extensive Monte Carlo 
simulation of the RFIM up to now by sampling 1280 disorder configurations for each 
lattice size and temperature and obtained via finite size scaling (for h/T = 0.35) 



(25) 



These are values for the binary hi = ±h distribution of the random fields, however, 
those obtained for a continuous (Gaussian) distribution are not significantly different 



V = 


= 0.56 


±0.03 


V = 


= 1.00 


±0.06 


V - 


= 1.6 


±0.3 


1 = 


= 2.3 


±0.3 


P - 


= 0.00 


±0.05 


a = 


= -1.0 


±0.3 



233|| . In addition Dayan et al. ^7j performed a real space renormalization group 
calculation that gave identical values for 77 and Tj as those listed in (^). Finally 
an extensive high-temperature series expansion by Gofman et al. ||105|| find in three 



dimensions 7 = 2.1 ± 0.2, concurring with the Monte Carlo data (p5[) and 7 



z/(4 — fj) = 27, the last equality being a consequence of (p^), which they showed 
to hold also in 4 and 5 dimensions. Moreover, they demonstrated an even stronger 
relation for the amplitude ratio to hold: 

A = hm -j#^ = 1 , (26) 

where Xdis = -Z^'^[(5'i)^]av is the disconnected susceptibility. The Monte Carlo data 



can be analyzed in a similar way ||233|| and agree with (^). Thus there are strong 
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indications for equation p^ ) to be correct and hence for the conjecture of two 
exponent scahng to be true. 

A closer inspection of the hst of exponents (pS]) shows some pecuharities: the or- 
der parameter exponent 13 turns out to be zero, which means that the magnetization 
jumps disco ntinuously to a non vanishing value at the critical temperature. This 
hints at a first order phase transition, a possibility that has already been suggested 
by Young and Nauenberg ||3U1| (see however the interpretation of X-ray scattering 
studies by Hill et al. ||122[ |). Note that /5 = is also found for a continuous dis- 
tribution ||233|| , which is important since in case of a binary distribution for larger 
field strength a tricritical point is predicted in mean field theory |jl|, which sepa- 
rates a first order (high fields) from a second order (low fields) transition line in the 
h — T-diagram. Usually, at a first order phase transition, one would expect phase 
coexistence at T^, which manifests itself in a typical multi peak structure in the 
probability distribution P{m) for the order parameter I^O], ^. In the Monte Carlo 
data no indication for such a scenario can be found |p33|| , instead P{\m\) shows a 
significant peak at a nonzero value for \m\ already at the transition T = T^ with 
P(0) — s> in the thermodynamic limit. 

Furthermore, at a first order phase transition the specific heat usually diverges 
with the volume of the system ||193|| , but the exponent a is negative, which means 
that the specific heat does not diverge at T^. By means of birefringence techniques 
Belanger et al. [0, |T3|, |216|| concluded from their experiments on Feo.47Zno.53F2 that 
a = (i.e. a logarithmic divergence of the specific heat). Only recently it was shown 
| 142|] that the same kind of experiments on Feo.85Mgo.i5Br2 are better compatible 
with a cusp like singularity of the specific heat and a = — 1, concurring with the 
value reported by in [p32|] see (|25|) . However, this value for a together with the other 
estimates in (^) would violate the modified hyperscaling relation (^). Moreover, 
Schwartz ||256|| derived an exact inequality 



2- a <iy{d -2 + 7]) 



(27) 



and accepting rj ^ 0.5 (since this result has a much smaller errorbar than u) it 
would imply u > 2, which is a rather large value. Indeed in a recent Migdal- 



Kadanoff study Cao and Machta [Q find such a large exponent u = 2.25, and they 
also report a = —1.37 and /5 = 0.02, consistent with (^). However, a value for a 
that is negative and large in modulus causes serious difficulties with respect to the 



Rushbroke-equality a + 2/3 + 7 = 2. Hence the random field enigma [ p61|| is still far 
from being resolved. 
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In analogy to the spin glass research activity (see previous sections) there has 
been considerable interest in the non-equilibrium dynamics also in random field sys- 
tems. It has been observed that diluted antiferromagnets fall into a metastable 
domain state if cooled in an external field B below the critical temperature Tc{B) 
(see the article of Kleemann for a review [13^1). This domain state has a finite corre- 
lation length that does not seem to increase with time, in contrast to the continuous 
aging phenomena observed in spin glasses (see section 3). Villain ||29(J|| and Grin- 
stein and Fernandez |l(J4| have predicted a logarithmic growth of the domain size 
(similar to what Fisher and Huse later predicted for spin glasses [0 ) , and numerical 
investigations of the domain growth in random field systems (e.g. [|223| , p21| , |212|| ) 
seem to be compatible with this prediction (see also the review of Chowdhury and 
Biswal [^). However, Nattermann and Vilfan [|191|| pointed out that DAFF map 
onto a RFIM plus random bonds, and the latter produce an enormous pinning force 
so that domain growth will only be observed on time scale up to 10^^ years, which 
explains the experimental situation. Moreover, these time-persistent domains are 
fractal objects and have been analyzed in Monte Carlo simulations of DAFF by 
Nowak and Usadel ||202| , P03|| . Recently also hysteresis effects in the RFIM at zero 
temperature have been studied |P58| , |76[| . 

Nattermann and Vilfan also predict that after switching off the external field in 
a DAFF the magnetization concentrated in the domain walls will decay according 
to M(t) ~ (logt)~^/*, which is compatible with experiments ||161| , |162|| . However 
the results of Monte Carlo simulations fit better to an algebraic decay [ pOCI| , |201 
an enhanced power law M(t) ~ Moexp{-A(logt/t')^} [|0|, |05|. 



or 



Very frequently it can be observed that by field decreasing experiments with 
a DAFF a "stable" domain state already occurs for field strengths B larger than 



the value Be below which the system orders antiferromagnetically (see [[L52|] for an 
overview). Thus in between the paramagnetic and the ordered phase a region for this 
domain state has to be inserted in the corresponding phase diagram (see figure ^. 
This state has been identified with an intermediate spin glass phase, characterized by 



dynamical freezing and lack of long range antiferromagnetic order [p.89| , |T^. From 
an experimental point of view (as well as for Monte Carlo simulations ||202|| ) it 
is still not clear, whether this intermediate regime corresponds to an equilibrium 
phase. However, only recently Mezard and Young ||182|| looked at the A^-component 
version of the RFIM in three dimensions and found in the limit N ^ oo that replica 
symmetry is already broken at the ferromagnetic transition at temperature Tf (which 
is usually an indication for a spin glass phase). In a subsequent investigation Mezard 
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and Monasson | 184 were able to show that above Tf a glassy phase appears in a 
temperature interval Tf < T < Ti,, Tf, being the temperature above which the 
usual paramagnetic phase is entered. Interestingly the ferromagnetic correlation 



length turns out to be finite, reminiscent of a result by Guagnelli et al. ||113|| , who 
study numerically the solutions of the mean field equations of the RFIM in three 
dimensions, and also very similar to the experimental finding of finite domain sizes 
in the intermediate domain state or spin glass region discussed above. Thus the 
phase diagram of RFIM (or DAFF) seems to be much richer than expected and in 
particular random field systems have much more in common with spin glasses than 
anticipated up to now. 

6 Quantum spin glasses 

The models discussed so far are all describing classical systems for which quantum 
fluctuations can be neglected. In most cases this is correct, namely as long as 
Tc > since critical fluctuations at the transition occur at a frequency huj <^ fc^T 
that is proportional to the inverse relaxation time r and thus approaches zero for 
T ^ Tc due to critical slowing down. Hence any finite temperature will destroy 
quantum coherence and the system will behave calassically. Very recently however, 
spin glasses began to enter the quantum regime ||244| . 

The interesting theoretical question is: What are the effects of quantum me- 
chanics on the physics of strongly disordered systems at zero temperature, where no 
heat bath is present and hopping over energy barriers is replaced by tunneling them 
quantum- mechanically. The best known and most studied example in this respect is 



the zero temperature metal- insulator transition ||165|| . The renewed interest in spin 
glasses in the quantum regime was kindled by a series of recent experiments ||297| , |298| 
on the dipolar Ising magnet Li2:Hoi_2:YF4, where Tc was driven down to zero by the 
application of a transverse magnetic field F (see the phase diagram depicted in figure 
1^). Experimental realizations of quantum spin glasses are already known for more 



than 10 years: the so called proton glasses p7[ pi5|| , which are random mixtures 
of ferroelectric and antiferroelectric materials like Rbi_a;(NH4)2.H2P04. Here the 
proton position, describable by an Ising spin variable, tunnels between two energy 
minima with a fixed frequency modelled by a transverse field acting on the spins. 
In the transverse field Ising magnets mentioned above, however, it became possible 
to study the zero temperature phase transition occuring for critical transverse field 
strength Tc by simply tuning the external field strength. This quantum phase tran- 
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sition lies within a different universality class than the usually studied (classical) 
spin glass transition at finite temperatures and it turns out that their properties 
differ significantly. 

The above mentioned experiments can be described by the model Hamiltonian 
of an Ising spin glass in a transverse field ||297|, |298| 



H = -Y. J^^^M - r E ^r f-^ E ^'] ' (28) 

where the ctj are Pauli spin matrices, F is the strength of the transverse field and h 
is a longitudinal magnetic field used to define magnetic susceptibilities but usually 
set to zero. Otherwise this model is equivalent to the EA-model in d dimensions 
defined in (|I]). Obviously, for F = the quantum-mechanical Hamiltonian (^) is 
diagonal in the z-representaion of the spin operators, which in this case can simply 
be replaced by their eigenvalues ±1 (after rescaling the couplings) giving exactly 
the classical EA-model (|l|). In this way the transverse field introduces the quantum 
mechanics into the spin glass problem and the value of F tunes the strength of the 
quantum fiuctuations. At zero temperature and F = the system described by ( |28D 
will be in its uniquely determined ground state, which is identical to the classical 
ground state of (|l]). In this case one has {a^) = ±1 for all sites i and therefore 
Qea = [{(^i)'^]av = 1 5 where (■ ■ ■) means the quantum-mechanical expectation value. 

If we switch on the transverse field (F > 0) the Hamiltonian (^) is not diagonal 
in the 2;-representation any more and its ground state will be a superposition of 
the classical ground state plus various excited states, which describes the quantum- 
mechanical tunneling at zero temperature between the local energy minima of the 
classical Hamiltonian. Furthermore | (erf) | < 1 since the transverse field tries to 
align the spins in the x-direction and therefore qea = [{(^i)^]a.v < 1- Increasing F 
diminishes the EA-order parameter qea and for some critical value F^ it will be zero: 
Qea = for F > Fc. This is the zero temperature phase transition we are interested 
in and obviously we cannot expect that its critical properties have anything in com- 
mon with the finite temperature classical spin glass transition discussed in section 
2. 

In order to describe this zero-temperature transition one introduces a quantity 
measuring the distance from the critical transverse field strength (at T = 0) 5 = 
(F — Fc)/Fc. If one assumes a conventional second order phase transition one has 
Qea ~ I'^I'^ and xsG ~ I'^l"''^ fo^' S ^ 0. Spatial correlations decay on a characteristic 
length scale that diverges at the critical point as usual: C, ~ 1^1"*^ and these exponents 
defined so far would be sufficient to describe the static critical behavior of a classical 
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spin glass transition. However, at a zero temperature transition driven solely by 
quantum fluctuations static and dynamic quantities are linked in such way that the 
introduction of a characteristic time-scale (or inverse frequency) is necessary (see 
e.g. 0): 



^r-r-i^r 



(29) 



where z is the dynamical exponent. This will become more evident below, when we 
consider an equivalent classical model. 



Much work in the past has been devoted to the infinite range model [Q, |135| , |158| , 
287 , 49, 106, 289] and the phase-diagram in the F-T-plane, which looks similar to the 
one shown in figure |^: for low enough temperature T and field F one finds a transition 
line separating a paramagnetic phase from a spin glass phase. Recently Miller 
and Huse ||188|] and Ye, Sachdev and Read |^99|| focused on the zero-temperature 
critical behavior and calculated the critical exponents 7 = 1/2 (with multiplicative 
logarithmic corrections), (3 = 1 and zu = 1/2. Interestingly it seems that this 
quantum SK-model seems to fall into the same universality class as an infinite range 
metallic spin glass model incorporating itinerant electrons, which was investigated 
recently by Oppermann [pl4|| . Another proposition for a mean-field quantum spin 
glass that is exactly solvable was made by Nieuwenhuizen [|196|| via the introduction 
of a quantum description of spherical spins. 

In one dimension the situation is quite different. The transverse Ising chain 
with random bonds and/or random transverse fields can be mapped (see below) 
onto the McCoy- Wu model ||176| , |177| , |1 78|| , for which various analytical results are 
known p3|, |6g, |194|, 0, [18|, [IS^]. In particular D. Fisher [||] has shown within a 
renormalization group calculation that typical and average spatial correlations be- 
have differently: the typical correlation length diverges with an exponent v = 1 (as 
also found by Shankar and Murthy [ P6CI|| ), whereas the average correlation length 
diverges with an exponent v = 2, thus obeying the rigorous inequality 2/v < d = 1 
59| as an equality. Due to activated dynamics it turns out that ^^ ~ ex^{A^/^) and 



therfore z = oo. Furthermore it is predicted that /? = (3 — a/5)/2 '^ 0.38 and that, 
similar to the surface magnetization [|1 78|| , the bulk magnetization behaves nonana- 
lytically as a function of the field h already above the critical point due to Griffith 
singularities, giving rise to a divergence of the linear and nonlinear susceptibility 
already in the disordered phase. In a recent finite size scaling analysis of results 
obtained via Monte Carlo simulations and the use of the transfer-matrix formalism 
69|| discrepancies to this scenario were found. This may well be a consequence of 



the smallness of the lattice sizes studied (L < 16), since for this situation one might 
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expect typical rather than average results. 

In two and three dimensions, which are most relevant for the above mentioned 
experiments, no analytical results are known — apart from renormalization group 
calculations using the Migdal-Kadanoff approximation |3^, |6^ . For this reason ex- 
tensive Monte Carlo simulations have been performed recently in two dimensions by 
Rieger and Young [ZS^] and in three dimensions by Guo, Bhatt and Huse [ I15|| . Usu- 



ally the investigation of quantum systems via Monte Carlo methods are hampered 
by various deficiencies, the sign problem being the most notorious one in this respect 
(see [|168|| for a review). In studying the Ising spin glass in a transverse field however, 
one can exploit the fact that it can be mapped exactly onto a classical Ising model 
described by a real Hamiltonian. Using the Suzuki- Trotter formula [p82| 



one can 



easily show that the ground state energy of the d-dimensional quantum mechanical 
model (^) is equal to the free energy of a ((i + l)-dimensional classical model, where 
the extra dimension corresponds to imaginary time, i.e. 

-^^^^ = limr^o{,Tre-^^ 

(30) 

where the imaginary time direction has been divided into L^- time slices of width 
Ar {AtLt- = j3), and the effective classical action, S, is given by 

S = -Y.T.K^Mr)SAr) -T.llKS,{r)S,{r +l)(-Y.T.HSM)] , (31) 

T {ij) T i \ T i J 

where the Siir) = ±1 are classical Ising spins, the indices i and j run over the sites 
of the original d-dimensional lattice and r = 1, 2, . . . , L,- denotes a time slice. In 
equation (|3l|) , 

H = Arh (32) 

exp{-2K) = tanh(Arr) 

One has the same random interactions in each time slice. In order to fulfill the 

second equality in ( ^0|) precisely, one has to perform the limit Ar -^ 0, which implies 

Kij —>■ and K —>■ oo. However, the universal properties of the phase transition are 

expected to be independent of Ar so we take Ar = 1 and set the standard deviation 

of the Kij to equal K. Thus K, which physically sets the relative strength of the 

transverse field and exchange terms in (^81), is like an inverse "temperature" for the 

effective classical model in (|3l|) . 

One sees that the (d+l)-dimensional classical model ( pi]) should order at low 

"temperature" (or coupling constant K) like a spin glass in the d spatial dimensions 
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and ferromagnetically in the imaginary time direction. From this one concludes the 
existence of two different diverging length-scales in the classical model (0): one for 
the spatial (spin glass)-correlations, which is ^, and one for imaginary time (ferro- 
magnetic) correlations, which is ^r- Thus in the representation (|3TD the link between 
statics and dynamics in the original quantum model ( PSf ) becomes most obvious. 
Correspondingly, to analyze the critical properties of the extremely anisotropic clas- 
sical model ( PTD one has to take into account these two length scales via anisotropic 
finite size scaling [^ . 



Monte Carlo simulations of the classical model (^) are straightforward — it 
turns out that sample-to-sample fluctuations are significant, for which reason one 



has to do an extensive disorder average [|236| , |115|| . However, the finite size scaling 



analysis is complicated by the fact that due to the existence of two diverging length 
scales ^ and ^^ one has to deal with two independent scaling variables: as usual L/^ 
and in addition the shape (or aspect ratio) Lr/L^ of the system |3^. Thus, with the 
usual definition of a spin glass overlap Q = L^'^L^^ J2i,T 'S'"(T)5'f (r) for the classical 
system, the dimensionless combination of moments of the order-parameter gfav obeys 

g^K, L, U) = 0.5[3 - {Q')/{QYU ~ gU^L^^", Lr/L') ■ (33) 

In isotropic systems one has z = 1, which makes the aspect ratio constant to one for 
the choice L = L^- and in order to determine the critical coupling K^ one exploits 
the fact that gav{K, L, L) becomes independent of L ioi K = Kc (see e. g. |2^). In 



the present case of a very anisotropic system z is not known a priori and one has to 
vary three different system parameters to obtain an estimate for K^ and z (and other 
exponents). The following method ||236| , |1 15|| enhances the efficiency of such a search 



in a three-parameter space and also produces reliable estimates for the quantities of 
interest: In the limit L^ ^ L^ the classical {d + l)-dimensional classical system is 
quasi-one-dimensional, and in the limit Lj. <^ L^ the system is quasi-(i-dimensional 
and well above its transition "temperature" in d dimensions (which is even zero for 
d = 2). Therfore one has gav{x,y) -^ for y ^ and for y — > oo. Hence, for 
fixed X, ga_^{x,y) must have a maximum for some value y = yma.x{x). The value of 
this maximum decreases with increasing L in the disordered phase K < K^ (where 
5 = {Kc/ K — 1) > 0) and increases with increasing L in the ordered phase. This 
criterion can be used to estimate the critical coupling, as exemplified in figure (^. 
If one plots ga^{Kc,L,LT-) versus L^/L^ with the correct choice for the dynamical 
exponent z, one should obtain a data-collaps for all system sizes L. Finally one uses 
systems with fixed aspect ratio Lr/ L^ to determine critical exponents via the usual 
one-parameter finite size scaling. 
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Various scaling predictions can be made if one supposes a conventional second 
order phase transition to occur at some critical "temperature" Kf. for the classical 
model. Hyperscaling in this particular situation would imply [^ 



a 



v{d + z) . 



(34) 



As usual one has 7 = u{2 — rj), where r] is defined via the decay of correlations at 

criticality: 

C{r) = [{S,ir)S,+r{r)yU ~ r~(^+^-^+v) ^ 

(35) 
Git) = [{S,iT)S,it + T))U ~ r-('^+-2+W2. ^ 

and from (|3^) one gets via a + 2/3 + 7 = 2 the relation 2/5/z/ = d + z — 2 + 7]. The 
uniform magnetic susceptibility defined via 



Xf 



dh 



i^r 



-7/ 



(36) 



with respect to the quantum mechanical Hamiltonian (psD, is related to the inte- 
grated onsite correlation function of the classical model (^): xf ~ J2tG(t) and 
therefore 7/ = /3 — uz. Analogously the divergence of the magnetic nonlinear sus- 
ceptibility for the quantum system 



Xnl 



dh^ 



i^r 



(37) 



can be estimated via the spin glass susceptibility of the classical model | |236| , |115 
giving 7' = v{2 — f] + 2z). 

In the following table we list the results obtained so far in various dimensions 





d = 1[92] 


d= 1|| 


rf = 2||236| 


d = 3 115| 


rf = 3|3 


rf > 6 P991 |188| 


Z/ 


00 


~ 1.7 


1.50 ±0.05 


- 1.3 


- 1.4 


2 


V 


2 


~ 1 


1.0 ±0.1 


-0.8 


-0.87 


1/4 


V 


0.38 


~0.4 


-0.5 


-0.9 


— 


2 


7/ 


div. 


-2.3 


-0.5 


finite 


— 


finite 


7' 


div. 


— 


-4.5 


-3.5 


— 


0.5 



(38) 

The symbol div. in the second column means that first and higher derivatives of the 
magnetization diverge already in the disordered phase. The word finite means that 
in three dimensions and for d larger than the upper critical dimension the uniform 
susceptibility does not diverge at the critical point. The results in the second column 
|9^ were obtained within a renormalization group calculation, those in columns 3 



to 5 1^, p36| , |115|| with Monte Carlo simulations, column 6 shows the result of a 



Migdal-Kadanoff RNG calculation (in [^ also 7 has been calculated in this way 
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giving 7 = 1.2, which is far off the MC-value) and the last column depicts analytical 
results from mean-field theory [p99| , |188|| . 

Note that the results reported so far are for the zero-temperature critical be- 
havior of the quantum model (pH]). To make contact with the recent experiments 
mentioned above, which are done in the vicinity of the quantum critical point at 
finite temperature, one can exploit the following scaling prediction. For finite tem- 
peratures the scaling variable of e.g. the nonlinear susceptibility is ^^ ■ T since a 
finite temperature for the quantum system ( P^ ) implies a finite length L^ ~ T~^ in 
the imaginary time direction of the classical model (|3l|) . Hence one has 

xUT,6)^6-^'Xr,i{^), (39) 

with Xniix) -^ const, for x ^ and Xni{x) — ^ x~^ ^^^ for x ^ oo (in order to cancel 
the divergent prefactor at finite T if 5 — * 0). Thus one has for F = Fc 

Xni ~ T-T'/- (40) 

and analoguously e.g. xf ~ J'-^fIzv ^ Equation (|40|) implies in three dimensions a 
strong divergence of the nonlinear susceptibility (when approaching T = 0) with a 
power ~ 2.7 that is very close to the one for classical Ising spin glasses (~ 2.9 see 
section 2). This is in striking contradiction to the observation made in the experi- 



ments mentioned in the introduction |p98|| and in order to clarify this discrepancy 
further work is necessary. 

Another very important issue is the effect of Griffith singularities on the prop- 



erties of quantum spin glasses. As already observed by McCoy and Wu |176|, |178[ 



for the one-dimensional quantum model (corresponding to the 2-dimensional clas- 
sical model with layered disorder) their existence lead to divergencies in the mag- 
netic susceptibility already in the disordered phase near the critical point (see also 
1 2601 , |195| , p^). Thill and Huse [p86|| recently presented a droplet theory for quan- 



tum spin glasses, where they also predict algebraically decaying correlations in this 
region. This implies that the effect of Griffith singularities on the time-persistence 
of autocorrelations is much stronger for quantum systems at zero temperature than 
their effect for classical models at finite temperatures, where only an enhanced power 
law is predicted ||222|| . Preliminary results of Monte Garlo simulations of the two- 



dimensional transverse field Ising spin glass at zero temperature give already strong 



support for this scenario to be correct [239 



Finally we should mention that also the more general case of spin-| Heisenberg 
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model with quenched disorder, described by the Hamiltionian 

H = Y: Mcr:[cT] + <5,,(araJ + afa|)} - ^/^ " ^ (41) 

is of great interest, in particular because e.g. the copper-oxygen layers of high-Tc 
superconductors are supposed to be good physical realizations of such a model in 
two-dimensions ||245|| . The additional disorder-parameters 6ij introduce a random 



anisotropy without destroying the XY-symmetry and the limit 6ij -^ reproduces 
the Ising case considered above. For vanishing external field h = the strength of 
the disorder {[Jij]a.v — [Jij]'i-v) triggers the zero-temperature quantum phase transition 
of the model pT]). For the isotropic case 6ij = 1 this scenario has been studied in two 
dimensions by Sandvik and Vekic ||245|| very recently with Monte Carlo simulations. 
Unfortunately, due to the notorious sign problem occuring in (Quantum)-Monte 
Carlo simulations, only the non-frustrated case {Jij > 0) has been considered. On 
the other side, much progress could be made in the study of the XXZ-chain with 
quenched disorder (i.e. model (^) in one dimension [^2|, |116| , p42|| ), which is well 



suited for exact diagonalization studies. 

7 Other models and conclusion 

Althought we have confined ourselves to Ising spin glasses in this review, we finally 
would like to mention some of the other spin glass like systems currently under 
investgation. A concise review on this matter has already been given by Young, 
Reger and Binder |p02|| . 

Many experimental spin glasses that most convincingly show features of a spin 
glass transition at finite temperatures are Heisenberg-like systems (e.g. CuMn and 
Eua;Sri_^.S) described by the (classical) Hamiltonian 

H = J2■J^JS^S^^ (42) 

where S_^ are 3-component vector spins. There seems to be ample numerical evidence 
rr^] for references) that this model (with short range interactions) does not 



see 



possess a finite T spin glass transition — which is an unsatisfactory discrepancy with 
the experimental situation (in four dimensions indications for a finite-T transition 
were found in a recent Monte Carlo study [^]). However, it has been demonstrated 



by Matsubara et al. [|174|] via numerical simulations using a hybrid Monte Carlo spin 



dynamics method that an anisotropy term modelled by 



H.n^s. = E E DS'srsf (43) 

(ij) a^P 
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induces a finite temperature phase transition (see also the defect-wall renormalization- 



group study of Gingras ||102|| ), which gives further support to the claim that all 
three-dimensional spin glasses with anisotropy should have a phase transition in the 
universality class of the short range Ising spin glass model ||6|. 

Another generalization of the Ising models considered here are so called Potts 
glasses (see | l302 | for a review), where each spin can take on p values instead of only 
p = 2 in the Ising case. It seems that in three dimensions for p = 3 this model 
shows a spin glass transition only at zero temperature [2^]. On the mean-field 
level it has been shown ||149|1 that the p-state Potts glasses and spin glass models 
with p-spin interactions show very similar static and dynamic behavior for certain 
values of p. The question, whether the situation is similar in finite dimensions. 



has been discussed in 2301, but due to the fact that on a finite dimensional lattice 



the definition of a model with p-spin interactions is to some extend arbitrary, it 
is not resolved yet. Furthermore it has been observed ||149|| that the dynamical 
saddle-point equations of both mean-field spin glass models (the Potts model and 
the p-spin model) have the same mathematical structure as that arising in mode- 



coupling theory for the structural glass transition ||139|| . This latter similarity also 



occurs for orientational glasses, (for a review see |Q), where the spins are replaced 
by quadrupolar moments. 

Very similar to the XY spin glass is the so-called gauge glass defined by the 
Hamiltonian 

H = Y1 Jij cos{ipi - ipj - Aij) , (44) 

where (f G [0, 27r[ are phase- variables and Aij G [0, 27r[ is a quenched random vector 
potential. This model has been introduced by Shih et al. ||263|| to describe granular 
superconductors in an external field. It has been discussed that the gauge glass 
(0) might describe some of the low temperature physics of vortex-glasses (see |p5|| , 
chapter 7 for a review), which is the proposed |p3| "true" superconducting state of 



high-Tc materials. Of particular interest in this respect is the existence of a gauge 
glass transition in three dimensions. Monte Carlo simulations [|130| , p27| , |228|| and 
defect-wall energy analysis ||101|] of the gauge glass found indeed a phase transition 
in three dimensions at finite temperatures. Recently, however, it was demonstrated 



| 38[| that screening effects might destroy this transition. Moreover, it is not clear 
whether the gauge glass and the vortex glass really belong to the same universality 
class ([^, chapter 7). The latter has gained much interest for the analytically 
tractable case in 1-1-1 dimensions, which has been studied numerically via Monte 
Carlo simulations only very recently |jll|, |7^. Descrepancies to existing theories 
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rekindled the interest in this model that still seems not to be fully understood. 

Finally we would like to mention the numerical investigation of a Hubbard model 
of interacting bosons within a quenched random potential |^|. Similar to the case 



of the quantum spin glass discussed in the last section this quantum mechanical 
Hamiltonian can be mapped onto a real classical Hamiltonian in d+1 dimensions. 
Monte Carlo simulations of this model in two dimensions | 275 , 293|| showed the 



existence of a transition from a super conducting phase to an insulating Bose glass 
phase. 

In summary we have seen that the application and usefulness of Monte Carlo 
simulations of strongly disordered systems has spread from Ising spin glasses and 
random field systems to very diverse fields including structural glasses, supercon- 
ductors and "dirty bosons". However, also the most "simple" spin glass model in 
finite dimensions, the EA- model with Ising spins, gains again increasing interest. 
Apart from many fundamental questions still being unsolved, like the existence of a 
transition within an external field, recent years' research activities began on one side 
to focus on quantum effects in spin glasses and on the other side to shift from crit- 
ical behavior to non-equilibrium dynamics. In both cases the main stimulus came 
from experiments and Monte Carlo simulations play an important role in learning 
to understand the newly observed phenomena. Thus, to conclude, spin glasses are 
far from being dead — they are still aging. 
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Figure Captions 



1. Contours of xsG = 1 (obtained via Monte Carlo simulations) for two-, three- 
and four-dimensional spin glasses in an external field compared to series ex- 
pansions |272|. From ref. |110|. 



2. Waiting time experiments: Relaxation of the tliermo- remanent magnetization 
for different waiting times tw The procedure mentioned in the text is sketched 
in the inset: the sample (CdCr1.7Ino.3S4) is cooled from above Tg = 16.7K 
down to 0.6Tg = lOK in a 15 Oe field, kept in the field during t^, then the 



field is cut off. From ref. 292 . 



3. Top: Decay of the remanent magnetization of (Fe2^Ni(i_2.))75Pi6B6Al3 as func- 
tion of time t for various temperatures {Tg = 22.6K) in a log-log plot. Bottom: 
Exponent m as a function of temperature obtained from a fit of the data of 
the upper figure to an algebraic decay Ms{t) oc t^™(^) for T < 0.97Tg. From 



ref. p7 | 



4. The autocorrelation function C(t, t^,) for the three-dimensional EA model with 
Gaussian couplings at T = 0.5 {Tg ^ 0.9 [^). The system size is L = 24 and 



the data are averaged over 128 samples. From ref. [|151|| , see [p34|| for similar 
data in case of the ±l-distribution. 

5. The correlation length S,{tw) of the two-dimensional EA-model with Gaus- 
sian couplings as a function of the waiting time t^ in a log-log plot. The 
straight lines are least-square fits toan algebraic growth law ([TB|) with expo- 
nents a{T) = 0.044, 0.066 and 0.081 for temperatures T = 0.2, 0.3 and 0.4, 



respectively. Data taken from ref. [|237 



6. Phase diagram of the diluted antiferromagnet in a uniform field: 

H = J2{ij) ^iCjSiSj — BY^i^iSi with 50% dilution (e = or 1 with probability 
1/2). PM means paramagnetic phase, AFM means antiferromagnetic long 
range order and SG is the intermediate spin glass phase (as explained in the 



text). From ref. |202 



7. Phase diagram of the diluted dipolar coupled Ising spin glass LiHoo.i67Yo.833F4 
in the F-T plane. Filled circle follow from the dynamical behavior of the 
linear susceptibility ||297||, open circles from measurements of the nonlinear 
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susceptibility. We are interested in the region close to T = 0, F = 1. From 
ref. p98 |. 



The dimensionless cumulant (7av for the (2+l)-dimensional classical model ( |3TD 
for three different values of the coupling constant versus the system size Lr 
in the imaginary time direction. The system size in the space direction is 
L = 4 (o), 6 (+), 8 (□), 12 (x) and 16 (A). Since the maximum of ^^avl-^r) is 
roughly independent of L at K^^ ^ 3.30 one concludes that the latter value 
is the critical coupling constant. From ref. ||236 |. . 
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